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ABSTRACT 

We examine the Kepler light curves of V1504 Cyg and V344 Lyr, encompassing ~736 d at 1 min cadence. 
During this span each system exhibited ^^64 — 65 outbursts, including six superoutbursts. We find that, in 
both systems, the normal outbursts lying between two superoutbursts increase in duration over time by a factor 
~ 1.2— 1.9, and then reset to a small value after the following superoutburst. In both systems the trend of 
quiescent intervals between normal outbursts is to increase to a local maximum about half way through the 
supercycle — the interval from one superoutburst to the next — and then to decrease back to a small value 
by the time of the next superoutburst. This is inconsistent with Osaki's thermal-tidal model, which predicts a 
monotonic increase in the quiescent intervals between normal outbursts during a supercycle. Also, most of the 
normal outbursts have an asymmetric, fast-rise/slower-decline shape, which would be consistent with outbursts 
triggered at large radii. The exponential rate of decay of the plateau phase of the superoutbursts is 8 d mag~^ 
for V1504 Cyg and 12 d mag~^ for V344 Lyr. This time scale gives a direct measure of the viscous time 
scale in the outer accretion disk given the expectation that the entire disk is in the hot, viscous state during 
superoutburst. The resulting constraint on the Shakura-Sunyaev parameter, Ohot — 0.1, is consistent with the 
value inferred from the fast dwarf nova decays. By looking at the slow decay rate for superoutbursts, which 
occur in systems below the period gap, in combination with the slow decay rate in one long outburst above the 
period gap (in U Gem), we infer a steep dependence of the decay rate on orbital period for long outbursts. We 
argue that this relation implies a steep dependence of Ofcoid on orbital period, which may be consistent with 
recent findings of Patterson, and is consistent with tidal torquing as being the dominant angular momentum 
transport mechanism in quiescent disks in interacting binary systems. 

Subject headings: accretion, accretion disks - binaries: close - novae, cataclysmic variables - stars: individual 
(V1504 Cygni, V344 Lyrae) 



1. INTRODUCTION 

Cataclysmic variables (CVs) are semi-detached interacting 
binaries in which a Roche lobe filUng K or M secondary trans- 
fers matter to a white dwarf (WD). CVs evolve to shorter or- 
bital periods and show a "gap" between = 2 and 3 (where 
Ph = Porbitai/1 lu") during which time the secondary star 
loses contact with its Roche lobe and mass transfer ceases. 
Thus the binary becomes fully detached. At P/j = 2 the 
secondary comes back into contact with its Roche lobe and 
mass transfer resumes. For < 2 angular momentum loss 
from the binary is thought to be due solely to gravitational 
radiation. The CV subclass of dwarf novae (DNe) are char- 
acterized by their semi-periodic outbursts. SU UMa stars are 
DNe lying below the period gap that exhibit short, normal 
outbursts (NOs) and superoutbursts (SOs). We refer the time 
from one SO to the next as the supercycle. SOs show su- 
perhumps which are modulations in the light curve at periods 
slightly exceeding the orbital period. There are two further 
subdivisions within the SU UMa grouping: (i) the VW Hyi 
stars at long orbital periods, near P^ = 2, for which the decay 
rate is fairly constant during a SO, and (ii) the WZ Sge stars at 
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short orbital periods, a little greater than Ph = 1, which have 
less frequent, larger amplitude SOs, for which the decay rate 
decreases during a SO. DNe outbursts are thought to be due to 
a limit cycle accretion disk instability (Lasota 2001) in which 
material is accumulated in quiescence and then dumped onto 
the WD during outburst. During short outbursts in longer pe- 
riod DNe, a few percent of the stored mass is accreted, and 
during long outbursts a significant fraction ~0.2 of the stored 
mass is accreted. For the SU UMa stars, a SO is thought to 
accrete > 0.7 — 0.8 of the stored mass. Although the accretion 
disk is never in steady state during the limit cycle, it is close to 
steady state during SO, with the local rate of accretion inside 
the disk \ M (r) | decreasing linearly from a maximum value at 
the disk inner edge to zero at the outer edge. The accretion 
disk modeling has traditionally been done within the Shakura 
& Sunyaev (1973, hereafter SS) formalism, using two values 
for the a viscosity parameter, ahot for gas in the hot, ionized 
disk, and acoid for gas in the quiescent disk. 

There are two bright SU UMa stars in the Kepler field ex- 
hibiting a variety of temporal behavior that make them worthy 
of a detailed statistical study of their outbursts, VI 504 Cyg 
{Kepler ID 7446357; Ph = 1.67) and V344 Lyr {Kepler ID 
7659570; Ph = 2.1). These are members of the VW Hyi sub- 
division. To date the two light curves have amassed 736.4 d 
at 1— min cadence. Excluding gaps and bad data points, the 
light curves contain 1000431 and 1000345 data entries, re- 
spectively. Previous studies of the Kepler data on SU UMa 
stars have found quiescent superhumps in V344 Lyr (Still et 
al. 2010), presented numerical models of the long term light 
curve of V344 Lyr (Cannizzo et al. 2010; hereafter CIO), and 
studied superhumps, both positive and negative, in the long 
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Fig. 1. — Kepler light curve for V1504 Cyg at 1 min cadence, spanning ~736 d. The light curve covers six superoutbursts and 59 normal outbursts. The 
vertical red lines indicate the local maxima for the normal outbursts in which coverage permits a reliable determination, and with well-sampled decays down to 
~2 mag below maximum. 



term V344 Lyr light curve (Wood et al. 201 1). 

Statistical studies of DNe have been useful in delineating 
the long-term properties of outbursts, e.g.., the outburst du- 
ration, recurrence times, and quiescent intervals, and placing 
physical constraints on models of outbursts (e.g., Campbell 
& Shapley 1940, Sterne, Campbell & Shapley 1940, Bath & 
van Paradijs 1983, van Paradijs 1983, Szkody & Mattel 1984, 
Cannizzo & Mattel 1992, 1998, Ak, et al. 2002, Simon 2004). 
Several interesting studies of the SU UMa stars have been car- 



ried out. For instance, van Paradijs (1983) studied the varia- 
tion in outburst duration ti, with orbital period for 15 DNe 
spanning the 2 — 3 hr CV period gap. He found that short 
outburst durations increase with orbital period, whereas long 
outburst durations are relatively constant with orbital period. 
Therefore for the SU UMa systems, which lie at short orbital 
period, the ratio tt, {SO) /tt, (NO) is large. The relation of su- 
peroutbursts to normal outbursts for DNe below the period 
gap and the relation of long outbursts to normal outbursts for 
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Fig. 2. — Kepler light curve for V344 Lyr at 1 min cadence, spanning ^736 d. The light curve encompasses six superoutbursts and 58 normal outbursts. The 
conventions are the same as in Figure 1 . 



DNe above the period gap are part of a general trend; super- 
outbursts are just long outbursts in short orbital period DNe. 
This finding was amplified by Ak et al. (2002) using a larger 
sample. In this work we seek to extend previous studies by an- 
alyzing the high fidelity Kepler light curves of two SU UMa 
stars to investigate the properties of the outbursts. In Section 
2 we examine the outburst properties of the NOs and SOs, in 
Section 3 we discuss the results, in particular the scaling of 
the SO decay rate with orbital period, and in Section 4 we 
summarize our results. 



2. OUTBURST PROPERTIES 

2.1. Normal Outbursts 

The relatively small number of total outbursts observed in 
the two systems to date does not yet allow one to amass de- 
tailed frequency histogram distributions, as was done for in- 
stance by Cannizzo & Mattei (1992) for the ^700 outbursts 
observed up to that time in SS Cyg. However, the variations of 
the outburst duration ti, and intervening quiescent interval tq 
within a supercycle contain information which can potentially 
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guide theoretical models. For instance, in the thermal-tidal 
model for the supercycle in the SU UMa systems (Ichikawa & 
Osaki 1992, 1994; Ichikawa, Hirose, & Osaki 1993), one sees 
a monotonically increasing sequence of NO quiescent inter- 
vals leading up to a SO. Now with the exquisite Kepler data, 
this can be directly tested. The NOs of the fainter SU UMa 
stars such as V344 Lyr have been poorly characterized prior 
to Kepler due to their intrinsic faintness (see, e.g.. Fig. 1 of 
Kato et al. 2002). 

Figures 1 and 2 present the currently available long term 
light curves for VI 504 Cyg and V344 Lyr. To obtain the 
apparent Kepler magnitude from the original data, electron 
flux per cadence, a correspondencaj of Kp = 12 to 10^ 
cadence"^ was adopted. Each light curve encompasses six 
SOs; V1504 Cyg had 59 NOs while V344 Lyr had 58. Two 
systematic effects are worth noting: (i) There are gaps in the 
data due to monthly data downloads, quarterly 90° spacecraft 
rolls, and spacecraft safe modes, (ii) In addition, there is 
no long-term absolute flux calibration for Kepler (Haas et al. 
2010, Koch et al. 2010): (1) Target masks do not contain all 
the flux from the target. The aperture is chosen to maximize 
a varying S/N over the whole quarter. As Kepler orbits the 
sun, however, targets drift across the mask as a result of dif- 
ferential velocity aberration. Varying amounts of target flux 
fall outside of the pixel aperture over the quarter, so the tar- 
get flux tends to vary smoothly over the quarter due to this 
effect. Another contributing factor to this effect is focus vari- 
ation due to the thermal behavior following spacecraft slews, 
especially immediately after Earth pointings for data down- 
load, or safe modes. (2) After each Kepler rotation, the target 
falls on a different CCD chip with a different PSF and the tar- 
get is centered differently, e.g., perhaps on a pixel boundary 
rather than a pixel center. The 4 arcsec pixels undersample 
the PSF. This means that the pixel mask changes shape and 
size each quarter, with different amount of light losses. (3) 
With different shaped masks, the number and amount of con- 
taminating sources in the mask will differ across the quarter 
gap- 

The lack of long term flux calibration is evident in the ele- 
vated flux level seen in V1504 Cyg during 100 < t{A) < 200 
(i.e, Q2 of observations). In addition the secular trending of 
the flux calibration within a given quarter may account for the 
apparent ^50 d e— folding decrease of the quiescent level fol- 
lowing the first two SOs in V344 Lyr, both of which occurred 
just after a roll maneuver In view of these considerations, 
we limit our attention in this study to timescales rather than 
amplitudes. 

In a standard light curve plot of magnitude versus time, a 
straight line decay means that the flux / = /oe~*/^, which 
we would therefore refer to as an exponential decay. A de- 
cay which deviates from exponential could be concave up- 
ward (slower-than-exponential), concave downward (faster- 
than-exponential), or some more complicated shape. The 
decays of the NOs are faster-than-exponential, which obvi- 
ates theoretical studies that attempt to introduce a weak radial 
power law dependency into the SS a parameter in order to get 
precisely exponential decays (e.g., Ichikawa, Hirose, & Osaki 
1993; Cannizzo 1994). We may quantify the deviation from 
exponentiality in the NO decays by dividing the decay into 
bins, each spanning a small range AKp = 0.15 — 0.20, and 
calculating the locally defined decay rate tdid mag~^) within 
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Fig. 3. — The first normal outburst in V344 Lyr, with eleven horizontal line 
segments along the decay at equally spaced intervals of AKp = 0.15. The 
segments define the edge of the ten bins used in calculating the histogram 
of decays shown in Figure 4. The decay for each NO is divided into similar 
bins, using the local maximum of each NO as the starting point. The local 
decay rate within a bin is defined to be = {tj — tj ^i) / {Kpj — Kpj^i), 
where the index j denotes a bin edge, I < j < 10. 



each bin. Dividing the decay into a total of 10 bins for each 
star, we take AKp = 0.15 for V1504 Cyg and AKp = 0.20 
for V344 Lyr For demonstration purposes, the binning dis- 
tribution for the decay of the first NO in V344 Lyr is shown 
in Figure 3. Within each bin we compute a frequency his- 
togram distribution of td values for all NO decays. The re- 
sults are shown in Figures 4 and 5. The small numbers in each 
panel indicate the magnitude range in Kp measured from the 
local maximum for each NO. Thus, the upper panels show 
the decay rates computed near maximum light, and as one 
progresses to lower panels one is sampling the fainter por- 
tions. For both systems one sees a shifting in the centroid 
of the frequency histogram distribution to lower td values as 
the decay proceeds, equivalent to saying that the decays are 
concave downward. In both systems, the parts of the decays 
fainter than about 1 mag below NO maximum asymptote to a 
centroid median of ~'0.6— 0.7 d mag^^. Thus the concave- 
downward aspect of the NO decays is due in large part to 
the rounded maxima; by the time the NO has faded to within 
AKp « 1 mag of quiescence, the rest of the decay is closely 
exponential. 

Bailey (1975) discovered that the decay rate for NOs is cor- 
related with orbital period (see his Figure 1 for eight well- 
studied systems), and derived the empirical relation td = 0.38 
d mag^^ Ph. Smak (1984) was the first to show that this re- 
lation can be used to infer ahot — 0.1. Subsequent work 
confirmed and expanded on this finding (Cannizzo 2001a). 
Warner (1995b) looked at a larger sample than Bailey and de- 
duced the fit td — 0.53 d mag~^ For the two systems of 
interest in this study, for which Ph — 1.67 (V1504 Cyg) and 
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Fig, 4. — A histogram of the fast decay rates for the 54 normal outbursts 
in V1504 Cyg with well-defined maxima and well-sampled decays, as shown 
in Figure 1. The numbers in each panel give the magnitude range in AKp 
down from maximum corresponding to the interval within which the local 
decay rate is measured, thus brighter portions of the decay are in the upper 
panels. 

2.1 (V344 Lyr), respectively, the predictions for the rates of 
decay would be 0.64 and 0.80 d mag^^ from the original Bai- 
ley relation, and 0.82 and 0.99 d mag^^ from the Warner fit. 
The presently determined values ~0.7 d mag~ ^ in V 1504 Cyg 
and ~0.6 d mag~^ in V344 Lyr for the faintest bins are closer 
to Bailey's original expression. If one averaged over the entire 
decay, the values would be larger, closer to Warner's expres- 
sion. The precision of the current determinations are limited 
by the small number of total NOs in both systems. With time, 
better statistics will allow better determinations, as has been 
done for instance with SS Cyg (Cannizzo & Mattel 1998, see 
their Figure 6). 

Figure 6 shows the outburst durations tf, for the two sys- 
tems. Based on the magnitude range for outbursts in the 
two systems, we use a cut of Kp = 16 in VI 504 Cyg and 
Kp = 17.9 in V344 Lyr to define outburst onset and termina- 
tion. There is an obvious trend of increasing outburst width 
tb between consecutive SOs: during a supercycle, ti, increases 
by ^20—70%. Figure 7 shows the asymmetry parameter f^, 
defined to be the ratio of the rise time to the total outburst 
duration, for the individual NOs depicted in Figures 1 and 2. 
Figure 8 shows the quiescent intervals tq between successive 
outbursts. In both systems there is a general trend of tq val- 
ues that increase to a local maximum about mid-way between 
successive SOs. For the last half of the V1504 Cyg light curve 
and the first half of the V344 Lyr light curve this effect is quite 
pronounced. In is interesting that while the values of fr and 
tq from one NO to the next are fairly continuous, even across 
an SO, the ti, values show a sharp discontinuity across an SO, 



Fig. 5. — A histogram of the fast decay rates for the 51 normal outbursts 
in V344 Lyr with well-defined maxima and well-sampled decays, as shown 
in Figure 2. Conventions are the same as in Figure 4. 

resetting to a much smaller value. 

At roughly the time of the third SO in V1504 Cyg, strong 
negative superhumps appear (Wood et al. 2011, in prepara- 
tion) and persist through the end of the light curve. The pres- 
ence of negative superhumps indicates that the mass is not 
accreting primarily at the rim of the disk. Although there is 
not an obvious change in fr at this time, it is perhaps note- 
worthy that tq increases. This may indicate a strong drop in 
the mass transfer rate from the secondary star, which would 
increase the recurrence time for outbursts. In the magnetic 
systems such as AM Her, which lack accretion disks, one has 
direct evidence for the instantaneous mass transfer rate from 
the secondary stars in the long term light curves, and one does 
often see intervals during which the mass transfer rate reduces 
to a low level (e.g., Hessman, Gansicke, & Mattel 2000). 

2.2. Superoutbursts 

Figure 9 shows a composite of the six SOs in VI 504 Cyg, 
positioned in time so that f = 4 d corresponds to the crossing 
of Kp =16 during SO onset. The rate of decay during SO is 
~8 d mag~^. Figure 10 shows a composite of the six SOs in 
V344 Lyr, positioned in time so that f = 3 d corresponds to 
the crossing of Kp = 18 during SO onset. The rate of decay 
during SO is ~12 d mag~^. Within each system the SOs are 
quite similar in overall appearance and duration. 

3. DISCUSSION 

A summary of the parameters associated with the outbursts 
in the two systems is given in Table 1 . The behavior of the 
sequences of outburst durations tb and quiescent intervals tq 
during a supercycle are quite different. The tb sequence shows 
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Fig. 6. — The values of the outburst durations tt, for the two systems, 
taking a cut line Kpcut = 16 for V1504 Cyg and Kpcut = 17.9 for V344 
Lyr. The vertical dotted lines (red) indicate the superoutbursts. The values 
of ti, for the .superoutbursts are tb(d) = 12.1, 12.0, 11.9, 12.0, and 12.0 for 
V1504 Cyg, and tt(d) = 16.7, 17.1, 16.7, 16.9, and 16.0 for V344 Lyr. Only 
five SOs are quoted for each system because of a data gap at the start of the 
sixth SO in both systems. 



TABLE 1 
Summary of Parameters 



Parameter 


V1504Cyg 


V344 Lyr 


Comment 


NO: ii,(d) 


1.1-2.9 


2.5-5.0 


NO duration 


SO: 4(d) 


12 ± 0.07 


16.7 ±0.4 


SO duration 


NO:/r 


0.19-0.41 


0.17-0.40 


NO rise/duration 




3.2-20.4 


2.5-21.5 


Quiescent duration 


NO: td(d mag^i) 


0.7 


0.6 


NO decay rate 


SO: taid mag-i) 


8.1 ± 0.4 


12 ±0.6 


SO decay rate 



an increase between SOs in both systems. This is consistent 
with a secular increase in the mass of the accretion disk. The 
tq values, by contrast, increase to a local maximum roughly 
half way between SOs before returning to their initial small 
value by the time the next SO begins. Thus the tq values are 
roughly continuous across SOs. In both systems the variation 
of tq values between SOs appears to be inconsistent with the 
prediction of the thermal-tidal instability model, which shows 
a monotonic increase in tq between SOs (Ichikawa, Hirose, 
& Osaki 1993, see their Figure 1; Osaki 2005, see his Figure 
3). This indicates that the triggering radius for NOs does not 
move uniformly outward with each successive NO, but rather 
attains a local maximum and then recedes. 

The asymmetry of the NOs is similar in the two systems, 
/,. ~ 0.2 — 0.4. According to theory, the degree of outburst 
asymmetry reveals whether the outburst is triggered at small 
or large radii in the disk (Smak 1984; Cannizzo, Wheeler, & 
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Fig. 7. — The values of the asymmetry parameter = trisc/ifc for the 
normal outbursts in the two systems. Only values for which the times of 
outburst maxima can be reliably determined are shown. The same cut line 
values are used as in Figure 6. 

Polidan 1986, hereafter CWP; Cannizzo 1998, see his Figure 
6 and associated discussion). Outbursts triggered at large radii 
tend to have fast, asymmetric rise times (fr < 0.5), while 
those triggered at small radii are more symmetric (fr ~ 0.5). 
Smak referred to these two types as Type A and Type B, re- 
spectively, whereas CWP used the terminology "outside-in" 
and "inside-out" outbursts. The low fr values in both systems 
studied in this work are consistent with the Type A/outside-in 
outbursts. It is interesting that the one outburst which seems 
to be certifiably symmetric and therefore presumably inside- 
out — outburst 1 of V344 Lyr — occurred after an unusually 
long quiescence interval, which may have led to a more filled 
disk at the end of quiescence (compare this outburst with that 
shown in Figure 5b of Smak 1984 and Figure 3 of CWP). 
During this part of the V344 Lyr light curve there were strong 
negative superhumps (Wood et al. 2011, see their Figure 11), 
indicating a tilted accretion disk. This would have the effect 
of spreading the accretion stream impact over a large range 
in radii, thereby lengthening the time to build up to the next 
outburst. 

The mean durations of SOs for VI 504 Cyg and V344 Lyr 
are 11.9 d and 16.8 d, and the e— folding rates of the (slow) 
decay during SO are 8.1±0.4 d mag^^ and 12.0±0.6 mag~^, 
respectively, where the errors indicate the standard deviations 
of the six measurements. Thus the SOs span ^1.5 e— foldings 
in V1504 Cyg, and --1.4 e-foldings in V344 Lyr Given that 
the SO accretion disks are close to steady state, the depletion 
fraction fd in the stored disk mass during SO for the two sys- 
tems, 

/d~l-e^*^(so)/t.(SO)^ 
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TABLE 2 

Slow Decay of Superoutbursts in SU UMa stars 



20 30 40 
outburst 



Fig. 8. — The values of the quiescence durations tq between consecutive 
outbursts for the two systems. As with Figure 6, the vertical dotted lines (red) 
indicate the positions of the superoutbursts. The same cut line values are used 
as in Figure 6. 

is ^0.78 and ^0.75, respectively. (These numbers probably 
over-estimate the depleted mass fraction because they encom- 
pass portions of the SO fainter than the e— folding plateaus.) 
These depletion fractions are larger than the ~ 0.2 inferred 
for long outbursts in DNe above the period gap (e.g., Can- 
nizzo 1993), and stem from the larger ahot/acoid ratios in the 
SU UMa systems. The ratio td(SO)/td(NO) ~ 20 is consis- 
tent with the previous interpretation of fast decays being due 
to the action of a cooling front within the disk, and slow de- 
cays being the result of gradual accretion onto the WD, and 
therefore governed by the viscous time scale in the outer disk. 
The viscous and thermal time scales can be written in terms 
of the local Keplerian rotation frequency and the local disk 
height to radius ratio h/r (Pringle 1981). The viscous time 
scale 

t ^ 2 



and the thermal time scale 



tth 



id 



(3) 



differ by a factor h/r, therefore h/r 1/20. 

CIO present an equation for the rate of decay of a SO, taken 
from Warner (1995a, 1995b), which is based on the viscous 
time in the outer accretion disk 



td 



t -nda'^'^ P^'^rn''' 



(4) 



where ahot,-! = ahot/0.1 and m\ = M\/Mq. This 
e— folding time scale should approximate the decay rate, ex- 
pressed in d mag^^. Hence in a large sample of DNe the 



System 


decay rate (d mag ^ 


) Ph 


Reference 


SWUMa 


10. 


1.36 


1 


WXCet 


10. 


1.40 


2 


FLTra 


7.7 


1.41 


3 


VI 028 Cyg 


8.3 


1.45 


4 


HO Del 


7.1 


1.50 


5 


GO Com 


6.7 


1.58 


6 


NSV 4838 


10. 


1.63 


7 


VWCrB 


10.3 


1.70 


8 


EG Aqr 


7.1 


1.83 


9 


V503 Cyg 


9.1 


1.86 


10 
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prediction is that the rate of decay during SO should vary 
as -P"'^^. Our analysis shows, however, that in going from 
V1504 Cyg with = 1.67 to V344 Lyr with P,, = 2.1, the 
SO decay rate increases from 8 d to 12 d. One does expect 
object-to-object variations in NO and SO properties which 
can be large compared to overall, secular trends with orbital 
period, but given the dependencies in Equation Q and the in- 
ference from a large number of DNe that cthot — 0.1, there is 
no way to account for such a steep variation, td oc Ph^- 

It is worth pausing for a moment and looking at how the 
SO slow decay rates for V1504 Cyg and V344 Lyr fit into the 
larger picture. Table 2 lists the slow decay rate of SOs for SU 
UMa stars taken from the literature, most of which are derived 
from VSNET data (see Kato et al. 2009, 2010). We only show 
data for the VW Hyi subclass for which the slow decay rate is 
relatively constant over the course of the SO. In Figure 1 1 we 
plot the SO decay rates for the systems shown in Table 2, as 
well as VI 504 Cyg, V344 Lyr, WZ Sge and U Gem. For WZ 
Sge we only plot the value for the initial, steep rate of decay, 
4 d mag^^. As noted in the Introduction, the SO decays in 
the WZ Sge stars exhibit a strong deviation from exponential- 
ity presumably due to the accretion of virtually all the stored 
mass. As the disk mass decreases, one expects a concomitant 
increase in the viscous time, leading to a pronounced concave 
upward SO decay. Therefore only the very initial SO decay 
rate gives a valid representation for the viscous time within 
the disk corresponding to the stored mass. As regards DNe 
above the period gap, virtually none of the long outbursts ever 
observed exhibited enough dynamic range in their overall de- 
cay that one could reliably extract the rate of decay. Only for 
the unusually long, 1985 October outburst of U Gem could 
this determination be made (Cannizzo, Gehrels, and Mattel 
2002), hence the single data point above the period gap in Fig- 
ure 1 1. If one takes the overall td — Ph trend between WZ Sge 
and U Gem as a true physical baseline, td oc Ph^, the trend 
between VI 504 Cyg and V344 Lyr in log t^ — log P^ space, 
td oc Ph^^ is roughly the same, although the normalization is 
offset slightly. Given the vertical scatter in Figure 1 1 for the 
other SU UMa systems, this coincidence is fortuitous. 

The failure of Equation Q to account for the difference 
between the SO slow decay rate in V1504 Cyg and V344 
Lyr indicates that a more sophisticated expression may be re- 
quired. One obvious deficiency is that there is no dependence 
on ttcoid which controls the available amount of fuel in qui- 
escence to power a SO. The physics behind the viscosity in 
quiescence is unknown at present, and one could in principle 
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Fig. 9. — A composite of six superoutbursts in V1504 Cyg, positioned 
in time so that t = 4 d con'esponds to the crossing of Kp = 16 during 
superoutburst onset. The line segments (red) indicate the linear least squares 
fit, for the range of data indicated. The rate of decay for the six superoutbursts 
are as shown in each panel. The mean value is 8.1 d mag~^. 



have large variations in Ofcoid between systems at similar or- 
bital periods. Let us go back to the definition of the viscous 
time scale = {aiiot^)~^ (r / h)'^ . Since the controlling time 
scale in this situation is the slowest one, this expression is to 
be evaluated at the outer disk radius r^. Following Cannizzo 
(2001b), we consider the hot state of the disk by adopting a 
Rosseland mean opacity k = 2.8 x lO^'* g cm^^pT^'^-^ and 
mean molecular weight /i = 0.67. We can then write the vis- 
cous time as 



t„ 



5.7 d 



5/14 13/14 -8/7 



:,io 



•'hot, 



S(r<i) 



2000 g cm- 



-3/7 



(5) 



where mi = A/wd/IA/q, r^^io = ra/{lO^° cm), and E(rrf) 
is the surface density at the outer disk edge. 

Looking at this expression, we see that the Warner (1995a, 
1995b) equation akeady contained an implicit assumption 
about E, as it only contains ahot and Ph. To derive a more 
meaningful expression for td ~ in the context of SOs we 
need to take into account two additional pieces of information, 
namely that during a SO essentially all the mass accumulated 
during quiescence is accreted onto the WD, and thereby lost 
from the disk, and also that the accumulated mass depends on 
the critical surface density Smax, which depends in turn on 



Fig. 10. — A composite of six superoutbursts in V344 Lyr, positioned 
in time so that t = 3 d con'esponds to the crossing of Kp = 18 during 
superoutburst onset. The line segments (red) indicate the linear least squares 
fit, for the range of data indicated. The rate of decay for the six superoutbursts 
are as shown in each panel. The mean value is 12.0 d mag~^. 



acoid- This scaling is given by (Lasota 2001) 



13.4 g cm 



-0.38^1.14 



(6) 



Since the attainment of E^ax is the triggering mech- 
anism for a DN outburst, the quiescence disk mass 
^(7r/2)r^Einax(?'d, ttcoid) at the time of SO onset can be 
written as 



AM, 



cold 



4.2 X 10^1 gmr°-3V|}K-or- (7) 



Once the SO has started, the mass that was present in the 
non-steady state quiescent disk gets redistributed into a quasi- 
steady disk with A/(r) ~ constant, and therefore with a radial 
profile E(r) « r~^/^ (SS). Therefore the mass of the hot disk 
can be written as 



AA-fhot = / 27rrdrE(r, ahot)- 



(8) 



Since the ratio ahot/ctcoid is large in SU UMa systems, a sub- 
stantial amount of mass needs to drain onto the WD before 
the surface is depleted enough for the cooling wave to start. 

We now set E(r, ahot) equal to its standard functional 
form for a SS disk with gas pressure and free-free opacity. 



^hot 



Eq a 



-4/5 -3/4 



hot 



tered into Equation 



Equal 



This is the same form that en- 



ing the disk mass at the end 
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of the quiescent state with that at the beginning of the SO, 
i.e., AMcoid = AA/hot, gives 



(9) 



the relevant S value for Equation (|5]). Therefore we can 
rewrite Equation ^ as 



f - Sn 9 H Tn°-52 0.36 -1.14 0.44 

— »u.z a rrii acoid'^hot,-i'^<i,io- 



(10) 



This expression differs from Equation Q in that it contains 
«coid> reflecting the fact that the total amount of matter stored 
in quiescence enters into the viscous time scale in the resultant 
SO disk. 

In extreme mass ratio systems such as the SU UMa stars, 
the Roche lobe of the primary asymptotes to ^0.82 of the 
binary separation (Eggleton 1983). Also, the relation be- 

2/3 

tween orbital separation and orbital period is aio = 3.53P^ , 
where aio is the orbital separation in units of 10^° cm. If we 
further take the outer disk radius to lie at 75% of the pri- 

9 /3 

mary Roche lobe, we have rd^w = 2.19P,j . Since it is fairly 
certain that a in outburst is always close to 0.1, ahot been 
scaled to that value in the previous equations. On the other 
hand, the a value in quiescence is uncertain, and has been tra- 
ditionally adjusted by trial and error for each system. For SU 
UMa stars, workers have typically had to adopt much smaller 
acoid values than for DNe above the period gap. Given this 
uncertainty, let us take Ofcoid to be a power law in orbital pe- 
riod, ttcoid = ctcfiPjj^- Making these substitutions into Equa- 
tion ( [TO] l we obtain 

= 113 d mr'alfa^^\^\Pli-''+'-''^. (11) 

The dynamic variation associated with one magnitude (2.512) 
is 0.924 times that associated with one e— folding (2.718). 
Therefore if we associate the rate of SO decay, expressed in 
d mag~^, with the viscous time at the outer edge, which is an 
e— folding time, then Equation (111 needs to be corrected by 
a factor 0.924. 

By comparing the SO decay rates of V 1504 Cyg and V344 
Lyr, we infer a decay rate id = 0.924t^ = 12 d (P/i/2.1)i-^ = 
3.16 d P^'^- Furthermore, adopting mi — 0.6 for the SU 
UMa systems gives mi'^^ = 0.767, and therefore 



80.2 d f P, 



0.29+0.36^^ 3^ ^g^p^l.i 



(12) 



Solving for the constant and exponent in the acoid scaling 
gives Qico = 1-25 x 10^"* and = 4.2 — a steep depen- 
dence of acoid on orbital period. For V344 Lyr at Ph = 2.1 
we obtain Ofcoid — 2.8 x 10"'^, which accounts for why CIO, 
in their modeling of V344 Lyr, had to adopt the small value 
ttcoid = 0.0025. The inference ip = 4.2 is only valid, how- 
ever, if the trend between VI 504 Cyg and V344 Lyr were 
representative of a trend with P/j, which seems unlikely given 
the data for the other SU UMa systems. Based only on the SU 
UMa data, it would be more likely that ac,o has a large vari- 
ation between the two systems. As noted earlier however, if 
one takes the WZ Sge -;-> U Gem interpolation to be physically 
meaningful, then the VI 504 Cyg O V344 Lyr interpolation is 
fortuitously coincident. 

A steep dependence of acoid on P^ may have observational 
support from another line of argument. In a recent study of the 
aggregate properties of DNe, Patterson (2011, see his Figure 
1 1) finds that the recurrence time for SOs has a steep inverse 



relation with binary mass ratio ii.ccui(SO) cx g-2.63±o.i7^ 
where q = M2/M1. Since the secondaries in cataclysmic 
variables are, by definition, in contact with their Roche lobes, 
this implies that M2 scales with orbital period, and hence 



trocur(SO) OC P,^ 



2.63±0.17 



Given that (1) SOs accrete vir- 



tually all their stored mass AAfcoid, (2) the mass feeding rate 
into the outer disk from the secondary Mt is thought to be 
relatively constant at ^1Q~'^'^ Mq yr~^ for systems below the 
period gap (Kolb et al. 1998; Howell et al. 2001), and (3) the 
recurrence time (Cannizzo et al. 1988) 



'^rccur 

(SO) 



AM,, 



.Id 



\Mj 



OC r, 



,3.14 



a. 



-0.83 
cold ' 



(13) 



it is quite surprising that iiccur(SO) apparently has a strong 
inverse relation with orbital period. As one progresses to 
shorter orbital period, one expects more frequent SOs since 
the disks are smaller and the critical mass is reached more 
quickly. We are seemingly drawn to the same conclusion as 
previously, that Ofcoid must have a strong dependence on or- 
bital period. Howell et al. (2001; see their Figure 2, panel (a)) 
indicates a rough dependence d log | Af^ | log P/i ~ 1 for 
the SU UMa stars. If we again set Ofcoid P^, substitute Ph 
for rd in Equation ( 13 1, and solve for ?/;, we obtain iIj = 4.5. 

If the overall trend in td for long outbursts between WZ Sge 
and U Gem, which is roughly matched (coincidentally) by that 
between VI 504 Cyg and V344 Lyr, is real, then two indepen- 
dent lines of reasoning imply a very steep dependence of acoid 
on orbital period. Why should one have ahot — constant, 
while at the same time acoid P^, where -0 — 4.2 — 4.5? 
Current thinking on the physical mechanism for viscosity in 
accretion disks centers on the magnetorotational instability 
(MRI - Balbus & Hawley 1998), an instabihty based on the 
shearing amplification of a weak magnetic field. The non- 
linear saturation limit of the MRI in ionized gas is thought 
to explain why ahot — >^0.1 (McKinney & Narayan 2007a, 
2007b). However, as the disk cools to the quiescent state and 
the supply of free electrons drops precipitously due to the ex- 
ponential dependence of partial ionization on temperature be- 
low ~'10^ K, the magnetic field no longer couples effectively 
to the gas, and fails to provide an effective viscosity. Menou 
(2000) tried to model the complete accretion disk limit cycle 
model for DNe using only the MRI, and failed. He discussed 
the failure of the MRI in quiescence, and predicted that if tidal 
torques (a completely different physical mechanism from the 
MRI) dominate the angular momentum transport in the qui- 
escent accretion disks in the SU UMa stars, one would ex- 
pect an anticorrelation between iiccur(SO) and q. He presents 
data for six SU UMa stars showing an anticorrelation. (Pat- 
terson [2011] shows the anticorrelation in much greater de- 
tail.) Menou notes that theories of tidally induced spiral waves 
or shock waves indicate a strong correlation between effec- 
tive torque and q (Papaloizou & Pringle 1977; Goldreich & 
Tremaine 1980). Our finding acoid oc P^'^ seems to sup- 
port the view that tidal torquing is the dominant angular mo- 
mentum transport mechanism operating in quiescent accretion 
disks in interacting binaries. 

For completeness we note that previous workers have com- 
mented on the apparent necessity of small Ofcoid to obtain 
the large amplitude SOs and long recurrence times in the 
short orbital period DNe. For example Meyer-Hofmeister 
et al. (1998) found that the large required stored disk mass 
and outburst long recurrence time in WZ Sge mandated 
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ttcoid ~ 0.001 in the disk instability model. Meyer & Meyer- 
Hofmeister (1999) argued for a possible physical cause in the 
loss of coronal activity in the very cool atmospheres of the 
low mass K-type secondaries expected in such short orbital 
period systems, resulting in lower magnetic fields in the ma- 
terial entrained in the LI point as it feeds the outer accretion 
disk. However, in view of the findings of Menou (2000) show- 
ing the lack of MRl in quiescent disks, and the fact that if the 
MRl were operating the strength of the final field within the 
disk would be unrelated to the amplitude of the seed field, we 
find the tidal torquing arguments more persuasive. 

The acoid value strongly influences two observables, the 
decay time of the viscous plateaus, and the recurrence times 
for outburst. It is important to verify whether the steep acoid 
relation mandated by the WZ Sge — U Gem trend shown in 
Figure 1 1 is consistent with the observed recurrence times for 
the two systems. In other words, does the relation Ocoid = 
1.25 X lO-'^P^-^ applied to both WZ Sge and U Gem lead 
to reasonable recurrence times for both systems? For WZ 
Sge, Ph — 1.361 implies acoid ~ 4.56 x 10^"*. From 
Smak (1993) and Steeghs et al. (2007) we have mi = 1, 
rio = 1.1, and Mt = 2 x 10^^ g s~^. (Smak's value for 
Mt may be too low, as it is below the floor value expected 
solely from gravitational radiation [Howell et al. 2001].) Us- 
ing the value for the stored quiescent mass from Equation (7) 
gives AA/coid ~ 4.37 x 10^"* g. Thus the recurrence time for 
SOs from Equation (13) is ^ 66 yr, versus the observed recur- 
rence time for superoutbursts in WZ Sge (no normal outbursts 
are seen) of ^ 33 yr. For a standard value given by gravita- 
tional radiation induced mass transfer below the period gap 
Mt — lO~^*'Af0 yr~^, the recurrence time would be ^ 22 
yr For U Gem, Ph = 4.183 implies ckcoid = 0.051, From 
Smak (2001) we have mi — 1.07, m2 — 0.39, orbital separa- 
tion a = 3.53 X 10^° cm Pl'^{mi + ms)^/^ = 1.04 x 10" 
cm, and rio = 0.3aio ~ 3. From Equation (7), AAfcoid = 
1.52 X 10^'' g for the stored mass. Given a mass transfer rate 
typical of systems above the period gap ~ 2 x 1O~^M0 yr^^ 
we use Equation (13) to arrive at a recurrence time ~ 146 
d. The observed recurrence time for outbursts in U Gem is 
^ 120 d, and for long outbursts only it is ^ 240 d (since the 
outburst pattern basically alternates between long and short 
outbursts). The long outbursts are more relevant as regards the 
recurrence time scale obtained from Equation (13) since it is 
assumed most of the stored mass accretes. Thus for both WZ 
Sge and U Gem the theoretical and observed recurrence times 
for long outbursts agree to within a factor of ^ 2, indicating 
both systems can have the same scaling of ctcoid on binary or- 
bital period. Since the orbital periods of the two systems are 
quite different, their values of acoid ^"e quite different. 

The e— folding decay times td — 10 d for the mass of the ac- 
cretion disk during SO in VI 504 Cyg and V344 Lyr represent 
^ 10"^ orbital time scales in the outer disk = 27rr2~^ — 540 

s ^^^^lo^- In other words t^/t^ ~ lO'^. The latest MRl 
simulations of accretion disks, which should be applicable to 
the outburst state of DNe, are not yet sophisticated enough 
to be directly compared to SOs. For instance, Sorathia et al. 
(2011) perform global magnetohydrodynamic simulations of 
accretion disks which span a dynamic range of only a factor 
of four in disk radii. They also neglect the vertical depen- 
dence of the gravitational potential, which does not allow the 
effects of magnetic buoyancy in reducing the magnetic pres- 
sure within the body of the disk. In their simulations they find 
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Fig. 1 1. — A comparison of superoutburst slow decay rates of other SU 
UMa stars, from Table 2, with VI 504 Cyg and V344 Lyr. Also shown are WZ 
Sge (initial superoutburst decay rate only) and U Gem. The vertical dashed 
lines indicate the 2 — 3 hr orbital period gap for DNe. 

td/t^ ~ 10^ for the depletion of accretion disk mass. Never- 
theless their results are quite interesting in the context of our 
work because they show that the Shakura-Sunyaev a values 
determined from the local r and components of the mag- 
netic field are about an order of magnitude lower than those 
determined from the averaged viscous time scale in the outer 
disk (equivalent to the global decay rate of the disk mass), 
which is the more relevant determination in the context of 
comparison with DN decay rates. King, Pringle, & Livio 
(2007) noted a factor ^ 10 discrepancy between observation- 
ally determined ahot values and those calculated numerically 
from MRl simulations (the MRl values being too small). Al- 
though the local shearing box simulations cited by King et 
al. have been strongly criticized because of convergence is- 
sues (Fromang & Papaloizou 2007), the global calculations 
of Sorathia et al. also give ahot — 0.01. If the general trend 
between two methods of ahot determinations by Sorathia et 
al. are confirmed in more sophisticated global MRl calcula- 
tions with larger routcr/'^inncr and realistic vertical stratifica- 
tion, that may resolve the current apparent discrepancy be- 
tween theory and observation. 

4. CONCLUSION 

We have analyzed long term Kepler light curves of VI 504 
Cyg, containing 65 outbursts, and V344 Lyr, containing 64 
outbursts. Each system showed six superoutbursts. The find- 
ings are: 

(1) The NO decays are faster-than-exponential, with most 
of the deviation occurring near maximum light. Near quies- 
cence the decays are close to exponential, ^0.7 d mag^^ for 
V1504 Cyg and ^0.6 d mag~^ for V344 Lyr These rates are 
in line with that expected from the Bailey relation. Our values 
are based on frequency histogram distributions of only ^50 
NOs and should improve as more outbursts are accumulated. 

(2) The NO outburst durations tb increase by a factor 
^1.2—1.9 over the time interval between consecutive super- 
outbursts, resetting to a small value after a superoutburst. This 
is consistent with a monotonic increase of the accretion disk 
mass with each successive NO. 

(3) The quiescent intervals tq between normal outbursts 
show a general trend in both systems of increasing to a lo- 
cal maximum about mid-way between superoutbursts. (There 
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was also an anomalously long quiescence interval at the start 
of the V344 Lyr light curve, potentially the result of a tilted 
disk.) This behavior for tg is inconsistent with Osaki's ther- 
mal tidal instability model for SOs, in which one expects a 
monotonic increase in tq values accompanying the monotonic 
build-up in disk mass and angular momentum leading to a SO. 
If tq is correlated with the NO triggering radius ^trig, then 
rtrig moves outward with each successive NO only through 
roughly the first half of a supercycle, and then recedes. 

(4) The NOs in both systems are asymmetric, i.e., irise/ife < 
0.5, consistent with "Type A" outbursts (Smak 1984). 

(5) The inference of the steep relation Ofcoid cx Pj^ gives 
strength to the notion of tidal torquing as the dominant an- 
gular momentum transport mechanism operating in quiescent 
accretion disks in interacting binaries. 

If one considers the slow decay rate of the long, viscous 
outbursts spanning DNe from WZ Sge to U Gem, one infers 
an overall variation (x P^'^- It is unfortunate that only one 
long outburst in one DN longward of the 2 — 3 hr period gap 
exhibited enough dynamic variation in flux during its plateau 
phase that a decay rate can be reliably extracted for the slow 
decay portion. It is fortuitous, given the general scatter in 
DN outburst properties at a given orbital period, that the de- 
cay rate between VI 504 Cyg and V344 Lyr shows a similar 
law cx Ph^- These dependencies are much steeper than the 
po.25 expected from the previous theoretical expression, due 
to Warner (1995a, 1995b). By starting with the definition for 
the viscous time scale in the outer disk during superoutburst, 
and taking into account the scaling for material accumulated 
during quiescence, we derive a more physically motivated for- 
mula, and find that a steep dependence of the quiescent state 



value of a is required, acoid a Pt^- I" simple terms, the dra- 
matic increase in stored quiescent mass with decreasing or- 
bital period brings about a decrease in the viscous time scale 
of matter in the superoutburst disk, which is formed out of re- 
distributed gas from the quiescent state. Recent work by Pat- 
terson (2011) on the recurrence time scale in SU UMa stars 
for SOs also supports a steep dependence of Ofcoid on Ph. 
The long outbursts in systems above the period gap and su- 
peroutbursts in systems below the gap both represent viscous 
decays, with the absence of transition waves in the disk. It is 
somewhat counterintuitive that the recurrence time for these 
long, viscous outbursts varies as P^^'^, and yet their decay 
time varies as P^'^- The steepness of the aco\d{Ph) relation 
explains why accretion disk modelers have had to adopt such 
small Q!coid values for the SU UMa stars, compared to DNe 
above the period gap, and also probably accounts for the fact 
that the extreme mass ratio black hole transient systems un- 
dergoing limit cycle accretion disk outbursts, like A0620-00, 
have such long recurrence times for outbursts. 

The Kepler data have already provided important con- 
straints on the physics of accretion disks, in particular the 
Shakura-Sunyaev a parameter: CIO showed that an interpo- 
lation between ckcoid and ahot based on the degree of partial 
ionization of gas was needed to account for the shoulders in 
the SOs of V344 Lyr, and in this work we have shown that the 
scaling of SO decay rate with orbital period mandates a steep 
dependence of acoid on orbital period, enforcing the notion of 
tidal torquing as being the dominant viscosity mechanism in 
quiescent accretion disks. 

We acknowledge the contributions of the entire Kepler 
team. 
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